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1 Introduction 



Electroweak processes in the two-nucleon systems invite detailed studies for multiple rea- 
sons. From the nuclear physics point of view, these processes offer a valuable testing ground 
of the basic inputs of nuclear physics. In astrophysics, the precise knowledge of the pp fusion 
cross section is of crucial importance for building a reliable model for stellar evolution pi. 
Furthermore, experiments at the Sudbury Neutrino Observatory (SNO) (2|, |3|] to observe 
solar neutrinos with a heavy-water Cerenkov counter have made it extremely important to 
estimate the vd reaction cross sections with high precision. 

In this note we study fid capture: fi~ + d — > + n + n, using effective field theory 
(EFT). Our work is connected to the above-mentioned urgent need of accurate estimates of 
the vd cross sections, o v &. To expound this connection, we will first explain the standard 
nuclear physics approach (SNPA), see e.g. Ref.fi]]. This is a highly successful method for 
describing nuclear responses to electroweak probes. In this approach we consider one-body 
(IB) impulse approximation terms and two-body (2B) exchange-current terms acting on 
non-relativistic nuclear wave functions, with the exchange currents derived from a one- 
boson exchange model. The vertices in the relevant Feynman diagrams are obtained from 
a Lagrangian constructed to satisfy the low-energy theorems and current algebra ||], while 
the nuclear wave functions are generated by solving the A-body Schrodinger equation, 
H\^a) = E\^a), where the Hamiltonian H contains realistic phenomenological nuclear 
potentials. The most elaborate study of fid capture based on SNPA was carried out by 
Tatara et al. (TKK) || and by Adam et al. M. 

Now, the best available estimation of a U( i based on SNPA is due to Nakamura et 
al. (NSGK) §§, while that based on EFT is due to Butler, Chen, and Kong (BCK) g. 
Since EFT is a general framework (TO], it can give model-independent results, provided 
all the low-energy coefficients (LEC) in the effective Lagrangian, C e ff, are predetermined. 
C-eff considered by BCK, however, does contain one unknown parameter (L±a), which 
they adjusted to reproduce a u d obtained by NSGK. After this adjustment, the results of 
BCK are found to be in perfect agreement with those of NSGK. The fact that an ab initio 
calculation (modulo one free parameter) based on EFT reproduces o v & of NSGK extremely 
well offers strong support to the calculation based on SNPA. At the same time, it stresses 
the importance of carrying out an EFT calculation free from an adjustable parameter. An 
interesting possibility is to use fid capture data as input to control the unknown LEC. An 
immediate question, however, is whether this process is "gentle" enough to be amenable 
to EFT. The substantial energy transfer accompanying the disappearance of a muon can 
lead to a region of the final three-particle phase space in which the intrinsic state of the 
two neutrons receives such a large momentum that the applicability of EFT becomes a 
delicate issue. Let this unfavorable kinematical region be called the "dangerous" region. 
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The problem of the dangerous region is reminiscent of the difficulty one encounters in 



applying EFT to threshold pion production in N + N -> N + N + vr [y , [TJ, [TJ, gj . It will 
turn out (see below), however, that, unlike the pion production case, fid capture receives 
only a tiny fraction of contribution from the dangerous region, and therefore the theoretical 
uncertainty caused by the dangerous region is practically negligible. 

The calculation in Q used the power divergence scheme (PDS) pq[ . We employ 
here another form of EFT, the Cut-off Weinberg scheme (CWS). In CWS, the transition 
operators are derived from irreducible diagrams in heavy-baryon chiral perturbation theory 
(HBxPT), and the nuclear matrix elements are obtained by sandwiching them between 
initial and final nuclear wave functions obtained in SNPA. As discussed in [16], a next- 
to-next-to-next-to-leading order (N 3 LO) calculation based on CWS contains one unknown 
LEC, denoted by d R . Like L\a discussed by BCK, the parameter d R controls the strength 
of a short-range exchange-current term and, once d R is fixed from data, we can make a 
definite EFT prediction for o v &. We therefore investigate here the relation between d R and 
the \id capture rate, T^- Our study is essentially of exploratory nature, given the present 
limited accuracy (see below) of the experimental value of T^. 

Another important point concerning d R is that, as emphasized in |16[] , the strength 
of d R can be reliably related to the tritium /3-decay rate, T%. Thus, using the experimental 
value of Fg, which is known with high precision, one can determine d R and then proceed 
to make predictions on various two-nucleon weak-interaction processes, including the fid 
capture rate, pp fusion rate, and vd cross sections. We will present here the first estimate 
of the [id capture rate obtained in this approach; the pp fusion rate has already been 
discussed in [16], and the I'd cross sections will be reported elsewhere [17|. 



2 The capture rate 

Although fid capture can in principle occur from the two fid hyperfine states (S^ = 1/2 
and S^d = 3/2), the capture is known to take place practically uniquely from the hyperfine 
doublet state. Therefore, concentrating on this dominant capture, we refer to hyperfine- 
doublet fid capture simply as fid capture and denote the hyperfine-doublet fid capture rate 
by T^ d . The measured value of is rjjjjf = 409 ±40 s" 1 g|] and T^J = 470 ±29 s" 1 p|. 
We remark that a high-precision measurement of is being contemplated at PSI |2 



Denoting by L the orbital angular momentum of the two-neutron relative motion in 
the final state, we can write 



L=0,l,2.. 



where F^ d is the rate of fid capture leading to the L state. Here we shall be primarily 
concerned with T^J since it is this quantity that contains information about d R . The 
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contributions of T^ d (L > 1) are significant,^ 1 but their calculation is not expected to 



involve any major EFT-related issues. In general, due to the centrifugal force, the L > 
1 contributions cannot be too sensitive to short-range physics, which implies that chiral 
expansion for them should converge rapidly. Specifically, the L = 1 contributions are 
dominated by the axial-charge (AC) and El transitions, whose one-body operators (which 
are NLO in chiral counting) are well-known. The lowest order meson-exchange corrections 
(MEC) to the one-body operators come from soft one-pion-exchange (OPE), which is N 2 LO 
in chiral counting. These soft-OPE terms, dictated by chiral symmetry, are well known, 
and they are model-independent. For L > 2 states, within the accuracy of our evaluation 
only one-body contributions || have to be included. In our exploratory study, therefore, 
we concentrate on a detailed evaluation of T^J , and for T^ d (L > 1) we simply use the 
results obtained by Ref . || . 

Muon capture by the deuteron is effectively described by the current-current hamil- 
tonian of weak interactions 



H - Gv 



J d 3 iI a (f)J a (f) + h.c, (2) 
where the leptonic and the hadronic charged currents are 



L a (x) = i> u {x)^ a (l - 7 5 )i/> /Lt (x) and 

J«(x) = (V a - A a ) a=l (x) - i(V a - A a ) a=2 (x), (3) 

respectively, and Gy = 1.14939 x 10~ 5 GeV~ 2 |2l]]; a (a) is the Lorentz (isospin) index. In 
the center-of-mass system of the initial \x-d atom from which capture occurs, we can safely 
assume = Pd = 0. Consequently, the four-momentum transfer to the leptonic system, 
q a = (p u — p^) a , reads (q°, q) = {E u — m^, p v ). The /i-d capture amplitude is then given 
by 



G 



(f\H w \i) = ^_ d (0) l a C$ nn (-q,p; SlS2 )\j a (q)\H> d (s d )), (4) 
where the Fourier transformed currents are 



j a (q) = J d 6 xe-^ x J a (x), (5) 
l a = e lx ' 9 {v(Pu,s u )\L a (x)\n~(p^, Sfj)) = u v (p u ,s u )j a (l - Ts^Gv^)- (6) 

In Eq.([|), x I' M _ ( i(0) = 1/ \fna,Q is the IS wave function of the [i-d atom, where ao = (m^ + 
fnd) / m ii m d a i with a ~ 1/137.036 the fine structure constant.^ 2 ^d(sd) in Eq.(|j) represents 



#1 For example, Y^ d 1 w |r M£ j according to Ref.| 



#2 The small correction due to the finite size of the deuteron is taken into account in the actual calculation. 
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the deuteron wave function with the z-component of its spin Sd] ^ nn (—q,p; S1S2) represents 
the final nn wave function, with total nn momentum —q, relative nn momentum p = 
(pi — P2)/2, and the z-components of the neutron spins, si and S2- It is easy to obtain 

|Gy^(0)| 2 f d 3 p f d 3 p u 



fj,d 



4(2J^ + l)y (2tt)3 J (2tt)3 



E E 



2vr5(A£) 



1 



1) s d\Jfidi Sfj,d} 



(7) 



where = ^, AE 1 is the energy difference between the final and initial states, Vl^ = ^d(sd) 
and ^ nn = ^ nn (-q,p; sis 2 ). 



3 HB^PT Lagrangian and the hadronic currents 

To calculate the capture rate we adopt Weinberg's power counting rule. In HB^PT the 
leading order (LO) Lagrangian is given by 

Co = B[iv-D + 2ig A S-A]B-lj2 C A(Br A B) 2 + f^Tr(iAHA^ + &Tr( X +)(8) 

1 A 4 

with D„ = d„ + \ [£t, d^] - itfn^ - i^^ f , A M = \ {ft, d^} + - §£C M ft and 

x+ = tfx& + €xH> where = T ( V m + -4°) and C v = I T ( V m ~ ^) denote the external 
gauge fields. In the absence of the external scalar and pseudo-scalar fields % = m w> an d we 
define the pion field as £ = exp (i^^j. It is convenient to choose the four- velocity and 
the spin operator S^ 1 as = (1, 0) and S'' 4 = ^0, ^j. 

The next-to- leading-order (NLO) Lagrangian (including the "I/ttt-tv" terms) in the 
one-nucleon sector is given in |22j while that in the two-nucleon sector is given in |12 
Combining them, we can write the NLO Lagrangian relevant to our case as 

Cy^y v — nt 11 ' ( a 2 \ 

2m y + Cl Tr X+ + I 4c 2 - J (v ■ iA) 2 + 4c 3 iA ■ iA+ 

+ (2c 4 + -1—) [S», S u ] [iA MJ iA u ] - 1 [5", 51 f+) B 

\ 2m N J r 771AT / 

4idi 55 -ABBB + 2id 2 e abc e^ x5 v»A"> a BS X T b B BS 6 t c B + • • • , (9) 



where e i23 = 1, A M = ^A^, and /+, = - <9„£ M - i [£ M , f + £, ] (djl v - d u TZ^ - 

i [JZn, T^-i/] )£• We find it convenient to use the dimensionless low-energy constants c's and 
d's defined by 

ci.2,3,4 = ci 2.3,4, and rfi 2 = gA „ 9 rfi,2- (10) 

7T1AT 



#3 



Our definition of the pion field differs from that used in Ref. ]13] by a minus sign. 
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The values of these low energy constants, 6^2,3,4, are taken from Ref. p2| : 

c 1 = -0.60 ±0.13, c 2 = 1.67 ±0.09, c 3 = -3.66 ± 0.08, c 4 = 2.11 ±0.08 (11) 

and cq = Ky = 3.70. These values were determined at tree level (or NLO) in the one-nucleon 
sector, which correspond to N 3 LO in our two-nucleon calculation. 



3.1 The one-body currents 

The one-body currents can be obtained either by an explicit HB^PT calculation or 
by the Foldy Wouthuysen (FW) reduction of the well-known relativistic expressions. The 
former method requires one-loop diagrams which consist of vertices from Cq (for N 2 LO 
contributions) and those which contain one vertex from £1 (for N 3 LO contribut ions); also 
needed are the corresponding counter-terms from £2 and £3. We adopt here the FW 
reduction method for convenience. Since the range of t = q 2 = m fl (m^ — 2E u ) for \id capture 
is small (— m 2 <t< m 2 ), the t-dependences in the standard form factors, F^ 2 {t) and GA{t), 
give only less than 2 % effects, which can be reliably taken into account by expansion in t, 
F?(t) = 1 + \r\ + 0(t 2 ), F?{t) = k v + 0(t) and G A (t) = g A (l + \r\ + 0(t 2 )). Keeping 
the terms linear in t is consistent with HB^PT to the order we calculate. Some caution is 
required for the Gp(t) term, which contains the pion-pole contributions, 



G P (t) 
1m 7v 



2m N G A {t) 



mi, 



(12) 



where (5{t) is a slowly varying function. Comparison with the explicit HBxPT calculation 
up to N 2 LO H| leads to 



m = f«g*NN _ 1 2 m 2 +0 (<P\ 1 + [( _ 2 . „ L5) ± o. 3] % . 
g A m N 6 V A v / 



(13) 



The authors of Ref. || found that is reduced only by ~2 % when (3 increases by 10 %. 
Thus, limiting ourselves to the f3 = 1 case entails at most 0.4 % error. ^ 
The resulting one-body vector (IB) current components are 



V 1B (q) 



En 

i 



-iq-Ti 



-iq-Ti 



1 + -rt 

6 v 8m 



N 



4m 2 N 



Pi + §(1 + K V )q x 3j 
m N 



Am 2 N 



+ (1 + 2k v ) ( iBi x Vi - -q, ^ 



V 



;i4) 



* 4 If the deviation of /3 from 1 were important, we would have to include the one-body pseudoscalar term, 
PiB in Eq.(p^), as is necessary for the two-body P2B = P in Eq.([l8|). 
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where T i = \(rf ' — irf) and^j = (p-+pi)/2. The one-body axial-vector current components 



can be written for convenience as 



A<? B + 



mt 



(15) 



which defines Aib, where 

i 



-iq-rt 



iq-ri 



TUN 



<*i ■ Q 



8m 2 N 



°i [ 1 + ) + 



2 (Pi Si- Pi- ViPi) -^qvi-q + iqxpi 

4m 2 N 



16) 



The above equations correspond to N LO in HB^PT [24]. Apart from the mentioned t 



dependence of the form factors, the N 2 LO contribution is found to be negligible, ~ 0.1 %, 
indicating a rapid convergence. We therefore limit ourselves to N 2 LO for the one-body 
currents although, to be completely consistent, we should in principle include N 3 LO . 



3.2 Two-body exchange currents 

Since the evaluation of the two-body exchange current is the focus of this work we 
discuss the various parts of this exchange current in detail. We write the two-body vector 
and axial- vector currents as V^^i,^) and A% B (ki,k2), where fcj = p- — Pi is the mo- 
mentum transferred to the i-th nucleon. For the two-body vector-charge current we know 
^23(^1,^2) = 0(Q 4 ) pql , which therefore can be ignored. The spatial components of the 



vector current have a one-pion-exchange contribution of order Q [26]: 



V 2B (ki,k 2 ) 



-i(n x r 2 )- H 



g\ f^i (^2 • k 2 ) 



<?2 (<?i • h) 



mi 



k 2 



mi 



+ 



at ■ ki 



01 ■ k 2 



mi 



mi 



(k 2 - h) 



k\ 

0(Q A ) 



(17) 



where ~ 93 MeV is the pion decay constant, and we make use of the notation {j\ x r 2 ) 
(ri x t 2 ) x — i(r\ x t 2 ) v . Analogously to Eq.([l5|), we write the axial-vector current as p7], [U 



A a 



^2B + 



mi 



(18) 



with 
A° 2B (h,k 2 ) 



-\- ( 9A 

1 (n x T2 > ^Tf 2 



(Ti • k\ 



mi 



k 2 



2g A 
m N f% 



+(l<-2), 



C 2 + C 3 " ^ I T 



_k\a\- k\ 



1 2 
mi 



k\ 



(19) 
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P(h,k 2 ) 



9A 



L , 1 +• 



i _ -» ^ 1 _ 

-(fi x f 2 ) pi + 4c 3 r 2 fc 2 + (c 4 + -)(ti x t 2 ) (<?i x A; 2 ) 



(n x f 2 ) (a i x 



cx 2 • k 2 



rat 



2di(r! <?i + r 2 a 2 ) + d 2 {r x x f 2 ) (<?i X a 2 ) + (1 <-> 2) } , (20) 
<x 2 • k 2 



9Am\ 
2rn N f? 



8cit 2 



+ (!<-► 2) 



(21) 



Only one combination of the LEC, d\ and d 2 , is relevant for the fid capture process, 

12 1 

d R = di + 2d 2 + -c 3 + -c 4 + -. 

6 d o 



(22) 



Exactly the same combination of LEC's appears in triton /3-decay, pp- fusion and the solar 
hep process|ll|. Adopting the same strategy as in Ref. |i~6[j , we fix d B from T^(exp), the 
experimental value of the tritium /3-decay rate. 

To facilitate the calculations, we perform a Fourier transformation (FT) of the above 
two-body currents. To control short-range physics in performing FT, we introduce a Gaus- 
sian cut-off regulator 



S A (k 2 



exp 



2 A 2 



(23) 



where A is a cut-off parameter. It is to be emphasized that, although our calculation without 
regularization involves no infinities, we still need a regulator since EFT, by definition, is 
valid only up to a certain momentum scale. The regulated delta and Yukawa functions read 



yoA{m,r) 



d?q 
(2^)3 

d 3 q 



A 3 

o- exp 

(4vr)f 



A 2 r r 



q 1 + m 2 



(24) 



We remark that this is exactly the same regularization method as used in Ref. [O]. 

In performing FT, we need to specify the time components of the momentum trans- 
ferred to the nucleons. Energy conservation imposes the constraint: k® + k\ = — q° = 
m ii — E v . In our calculation we will adopt the so-called fixed-kinematics assumption (FKA) 
[11], where the energy transfer is assumed to be shared equally between the two nucleons, 
i.e. k® = k® = m>L 2 — , which naturally brings in the quantity = yjm 2 — (m^ — E v ) 2 /4. 
The uncertainty related to FKA becomes large as grows. The contribution from the 
large region, however, will turn out to be so tiny that the assumptions related to kf 
cause little uncertainty in our calculation. 
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4 The capture rate for the transition to the nn state 

The deuteron and the 1 5'o wave function may be written as 



ipd{r; s d ) 



4?rr 



/ \ ■ S 12 (f) 
Ud{r) + —j=-w d (r) 



Xi,s d £o,a, i>o(r) 



47rr 



uo(r)xa,o 6,-1 (25) 



with dr [urf(r) + w^(r)] = 1 and lim r 



sin 8a 



[cos pr + cot 5q sin pr] . Here 



Sizif) = 3c?i • r a 2 • r — d\ ■ oii X (0 is t ne Pauli spinor (isospinor), and 5q is the nn 
1 5q phase shift. To facilitate numerical work, we approximate AE as 



(26) 



where m = m n = 939.566 MeV is the neutron mass, M^ d = + m d = 1981.272 MeV. 
In our calculation we will neglect the O f^%jfr~) term since, as we shall show, the major 
contributions comes from the low p = \p\ region. Choosing the z-axis along p v , we write 
q = Pu = E v z. This simplifies the structure of the transition amplitudes as 



&o\f(q)\Vd(sd)) = 5 Sd ,oM u (Mel • M\*d(sd)) = 5s d ,\M x , 



(27) 



where e± = =f(£ ± w)/V^, &o = z, and A = ±1,0. We decompose the matrix elements into 
vector and axial vector current contributions, A4t,x = -A^aIX] — .M^aLA], and arrive at 



fid 



2tt 2 



dp2p 2 E 2 u [ 1 



E„ 



M, 



fid 



2M-i+Mo - Mi 



(28) 



Note that M-i = —(M+i[V] + M + i[A]), M Q = -M [A] and, to the order under consid- 
eration, Ait = — -Mt[A]. The matrix elements of the vector current are 



M\[V] = \V2 J^dr |< 
+ A(4x/2)f-|| 



^0 Ud Jo 



J2 W d \ fly 



s/2 J 2m N 



(!) , 
LO U\' [U d + 



W d \ . 2\iy - 1 



V2j Jl 4m2 



N 



1 

druQ I dx 

2 



x 32 w d 

JO u d 



xqrjf ( u d + ^ ) y L + J ( jf u d - ( v^'g + ^§= ) w d ) yf 



V2 

if 



L 2 r 

y - 3^ 



V2 



(29) 



where = jn(qrx) are the spherical Bessel functions, y„ = y n A ( ym 2 H ■^-q 2 -,r ) and 

yiA(m,r) = -r-§py A{m,r) and y 2 k{m,r) = -^r^-j^y^m, r). 
Using Eq.(ITH), we obtain for the axial current 



M t [A] 
M X [A] 



M t [A] 
M X [A] 



+ 



1- M <5a,o|^1 



a;M[i]-k1^o[i]+A^[^] , (30) 



9 



M t [A] 



M X [A] 



M[P] 
where 



V2g A 
V2g A 

■I 7T 

V2g A 
1 



dr \ —u§\u d - V2w d )ji - -^ruoiudjo + V2w d j 2 ) 

( - , * 9a\ , m ii -E v 
- c 2 + c 3 - -f- 1 



8m 2 N 



m N 



dru (ud - V2 w d )ji — , (31) 
r 



dr 



-, t 2 



^2 
P 

3m 2 N 



8m 2 N 



W d \ 



,(2) 



u d 



w d \ -A /TT 



(4V2) 



r 







— (1 - 5\fl)\q\u - (U d + -j=)]l 



^Voa ( c 3 + 2c 4 + ^ ) u ( u d jo - — |i2 



+- 



2/2A c 3 - c 4 



3 

.9-4 



"0 



V2 

w d jo - {u d 



(4V2) 



2m N f 2 Jo 



dr 



~ 2V1A , 

2cim w u {u d 

r 



V2w d )j! 



(32) 
(33) 



I 0\ 

= -Sx,o—u (u d - v2 w d )h 

o 

+ j^tio + 32 ) Mu' d - V2 w' d ) 



u (u d 



V2w d ) 



4^2 



(2j - ji)uow d , (34) 



j 2 A = (l-3<5 A , )i2, jn =jn(W), n«(r) = <(r)-^and u^(r) = ^( r )- 3 ^ + 3^. 
In the above expressions the curly brackets denote IB contributions, and for clarity we have 
suppressed the dependence on r in some equations. 



5 Results 

Table |l] shows T^J as a function of the cut-off parameter, A. As discussed, the short- 
range exchange current contribution depends on the single low-energy constant d R , see Eqs. 
d32|), (22), and d R determined from To(exp) is a function of A (see Ref.|16[|). We observe 
that the variation of T^J over the range of A under consideration is less than 0.7 s _1 . 
The d^-dependence in the table indicates the importance of the contribution of the short- 
distance exchange current. Without the d R term, T^J would change as much as 16 for 
A = 500 - 800 MeV. Thus, renormalizing the d^-term using T^(exp) reduces the variation 
of with respect to A by a factor ~ 20, leading to the practically A-independent behavior 
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A (MeV) 


d R 




500 
600 
800 


1.00 ±0.07 
1.78 ±0.08 
3.90 ±0.10 


254.7 - 9.85 d R + 0.159 (d R f = 245.0 ± 0.7 
261.1 - 9.09 d R + 0.132 (d R ) 2 = 245.3 ± 0.7 
271.0 - 6.76 d R + 0.070 (d R ) 2 = 245.7 ± 0.6 



Table 1: L = capture rate (in s 1 ) calculated as a function of the cutoff A. Also listed 
are the corresponding values of d R determined from r*j(exp) fl~6|| . 



of r^. Considering this stability we will hereafter only discuss the case corresponding to 
A = 600 MeV and d R = 1.78. 

The capture rate contains several interference terms, which are listed in Table || 
in a cumulative manner. We note that the axial charge (AC) plays only a minor role; its 



I> is' 1 } 




L = 




L > 1 


Total 




|GT| 2 


|GT±AC| 2 


|GT±AC±M1| 2 






1 1B I 2 


178 


177 


232 


138 


370 


|1B+2B| 2 


187 


186 


245 


141 


386 



Table 2: Cumulative contributions to (calculated for A = 600 MeV and d R = 1.78). 
The row labeled "IB" corresponds to the case that contains one-body contributions only, 
while the row labeled "1B±2B" to the case that includes both the one-body and MEC 
contributions. The three columns labeled "L=0" show contributions from the L = channel, 
with the contributions of the different transition operators displayed in a cumulative manner. 
The fifth column gives contribution from the L > 1 channels, as evaluated in TKK, Ref.Q, 
and the last column shows the sum of the L = and L > 1 contributions. 



destructive interference with GT decreases the capture rate by ~ 1 s _1 . Meanwhile, the Ml 
contribution interferes constructively with GT, increasing by ~ 59 Furthermore, 
the two-body MEC in the L = channel increases the capture rate by ~ 13 s _1 . 

Our final result for T^J = 245 s^ 1 in Table Q should be compared with TKK's 

lid 



result, r£j°(TKK) = 259 s~ l .W^ By adding the 1 < L < 5 contribution, r£f 1 = 141 s' 1 , 



fid 

calculated by TKK, we arrive at the total capture rate 

r^sse^- 1 , (35) 

#5 We have re-run the code of TKK using = 1.267. TKK's original result corresponding to qa = 1.262 
was r£r°(TKK)= 257 s _1 . 
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to be compared with TKK's result r^TKK) = (397 ~ 400) s -1 . 

As mentioned earlier, a primary question is whether ud capture process is "gentle 
enough" for applying HBxPT with reasonable confidence. As noted the "dangerous" region 
for HB^PT occurs when the two neutrons carry most of the final energy. To address 
this issue, it is useful to consider the differential capture rate, dT ^1 dE nn , where E nn = 
2(-\/ _|_ jp. _ m ^ i g tjjg energy of the final two-neutron relative motion. An equally 
informative quantity is the "cumulative" capture rate 

r„ d (E nn ) = / (dT^ d /dE' nn ) dE' nn . (36) 
Jo 

From these quantities we can assess to what extent \xd capture is free from the "dangerous" 
kinematic region. We show in Table ^ the matrix elements, A4+i[j4], .M+ifV], A4o[^4] and 



E 


M+i[A] 


M +1 [V] 


M [A] 


M t [A] 




0.0 


73.09 + 1.24 


14.68 + 0.53 


50.22 + 0.81 


0.79 - 0.23 





1.0 


20.88 + 0.38 


4.15 + 0.16 


14.26 + 0.25 


0.18 - 0.07 


91 


10.0 


2.59 + 0.12 


0.47 + 0.04 


1.82 + 0.08 


0.06 - 0.01 


231 


30.0 


0.49 + 0.05 


0.07 + 0.01 


0.39 + 0.04 


0.04 - 0.00 


244 


pmax 


0.056 - 0.003 





0.056 - 0.003 





245 



Table 3: Matrix elements calculated for representative values of E nn (MeV) and the cu- 
mulative L = capture rate for the case: A = 600 MeV and d R = 1.78. In each entry 
for the matrix element, the first number (preceding a "+" or "— " sign) gives the one-body 
contribution, while the second number gives the two-body contribution. 

.Mt[A], calculated for representative values of E nn , and for A = 600 MeV and d R = 1.78. 
Table || also gives ^ (^E nn ) . The graphical representation of r^ t( ^(£' Jin ) can be found in Fig. 
ffl. We learn from Table [|| that the matrix elements decrease quite fast as E nn increases, 
a feature that can be easily understood as follows. The 1 S$ nn radial wave function is 
proportional to = -j- [(pcot^o) 2 + p 2 ] 1 ■ Since the nn scattering length is very large, 
p cot 5q diminishes rapidly when the nn relative momentum p gets small. The examination 
of Table || also reveals that the one-body amplitudes decrease more quickly than the two- 
body amplitudes. This is a consequence of the softness of the deuteron wave function, 
which cannot supply high momentum transfers needed for producing large values of p. As 
a result, the contributions from high E nn - where the applicability of EFT is questionable 
- is negligible. For instance, the contribution to T^J from E nn > 30 MeV is just 1.1 s -1 , 
and that from E nn > 50 MeV is less than 0.1 

We now can make a rough estimate of the theoretical error associated with this 
calculation. Uncertainty related to the Gp term, Eq.([l2]) (or j3) is ~ 1 s , while uncertainty 
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400 




E nn (MeV) 

Figure 1: Cumulative fi-d capture rate (in s _1 ) calculated for A = 600 MeV and d R = 1.78. 
The dashed line gives the L = contribution, T^°(E nn ), while the solid line shows the 
total contribution, T^iEnn) = r^J°(-EVm) + r^p(-Enn)- The empty and solid circles for 
the values at E nn = £™ n ax ~ 102 MeV, for L = and L > 0, respectively. 

reflecting the A-dependence is less than 1 s -1 ; uncertainty in T^(exp) (or that in d R for a 
given A) can affect Y^d at the level of 1 s . If we assign a rather conservative error, 2 s , 
to the L > 1 contributions obtained in Ref. ||, the overall uncertainty in our estimate 
becomes 4 s _1 or 1 % in the total capture rate. 

As mentioned, there is a serious disagreement between the two measured values of T^. 
Our theoretical result is consistent with Y^iexp) in Ref. [lg| . In the present exploratory 
study we have not considered radiative corrections |28||, which are expected to be smaller 
than the existing uncertainty in Y^iexp). When the planned precision measurement of the 
T^d at PSI 0| is realized, the issue of radiative corrections should certainly be addressed. 
The EFT approach as described here will provide a useful tool for this purpose as well. 
Once the accuracy in r^(exp) is significantly improved, we will be able to use fid capture 
to determine the low energy constant d R , a quantity critically important for the accurate 
evaluation of the vd cross sections used in the analysis of the SNO experiments. At present 
the tritium /3-decay is a much more accurate source of information on d R than \xd capture, 
but it is hoped that in the near future T ^ will provide an independent constraint on d R . 
We consider this redundancy extremely important. 
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